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Abstract 

The volume operator plays a central role in both the kinematics and dynamics of canonical 
approaches to quantum gravity which are based on algebras of generalized Wilson loops. We 
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1 Introduction 


Researchers in gravitational physics these days can look back on a ten-year long effort of 
quantizing the theory canonically in terms of a set of “new” connection variables [1]. These 
developments have led to many new insights, but at the same time have not been free of 
Irrungen und Wirrungen, and the process is by no means finished yet. 

From a technical point of view, the problem of representing non-polynomial quantities 
in the quantum theory (a central difficulty in ADM-type quantization approaches) in this 
formulation has now been recast into that of diagonalizing the operator det E, where det E 
is up to a sign the classical determinant det g of the Riemannian metric g ab on spatial slices. 
The fact that non-polynomiality can be re-expressed as “polynomiality modulo powers of 
y / det g ” is not specific to the connection formulation (see, for example, [2]). However, what 
distinguishes this new approach is the fact that (certain functions of) det g have well-defined 
finite, self-adjoint operator analogues. 

This happens because in the quantization the classical conjugate variable pairs ( A l al E “) - 
like in Yang-Mills theory - are represented by multiplication and differentiation with respect 
to A l a . This implies that suitably quantized versions of the classical determinant 

det E(x) = ^r labc e i ^E^x)E b j (x)E c k ( x) (1.1) 

are third-order differential operators. Moreover, if the quantization is based on one-dimensio¬ 
nal flux line states, as is the case in the so-called “loop quantization” schemes [3,4], and also 
in lattice-discretized versions of canonical quantum gravity [5], R-fiux is quantized. This is 
ultimately responsible for the finiteness of quantities like the volume operator, which is the 
quantization of the integral d 3 x Vdet ~E. for a spatial region 1Z. 

The framework we will be using in the following is not the one originally proposed in 
[1], since the reality conditions that have to be imposed on its si(2, (D)-valued connections 
A l a (x) seem to be intractable in the quantum theory. Instead, we will use a closely related 
version in terms of real, sn(2)-valued connection forms [6], which avoids this difficulty. In 
this formulation, the Hamiltonian constraint function is non-polynomial. 

The subject of this paper is the introduction of a method for simplifying the diagona- 
lization of the volume operator (and potentially other operators relevant in loop quantum 
gravity), by exploiting symmetries of the Wilson loop states that form the Hilbert space it is 
defined on. Our discussion will take place in the discretized version of the theory on a cubic 
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lattice. The discrete volume operator one can define on the lattice [7,8] is closely related to 
the finite volume operators one obtains after regularization in the continuum theory [9-14]. 
The type of intersection of Wilson loop states that can occur on the lattice is at most six- 
valent, and therefore not the most general from the point of view of the continuum theory, 
but also in this case it may well be sufficiently generic from a physical standpoint. 

The operator detH occurs in a variety of contexts. It was originally conceived as an 
ingredient in the definition of the volume operator, i.e. the quantization of J .^ d 3 x \J det E 
mentioned above [9,10]. Defined on the kinematical Hilbert space of the loop representation, 
this kind of geometric operator enables one to associate well-known geometric properties with 
Wilson loop states, or possibly coarse-grained ensembles of such states [9]. 

With the return to real connection variables in canonical quantum gravity, knowledge of 
the spectrum of det E has become a necessary prerequisite in the study of the quantum Hamil¬ 
tonian constraint. The fact that the quantum dynamics can still be made well-defined, in spite 
of the non-polynomiality of the Hamiltonian, was first demonstrated in the context of the lat¬ 
tice theory [15]. A quantum Hamiltonian for the continuum theory was constructed in [16]. In 
order to get rid of the inverse powers of det E that occur naturally, an identity was employed 
by which the dreibeins e l a (which are non-polynomial functions of the inverse, densitized 
dreibeins Ef that constitute half of the basic canonical variables) can be re-expressed as the 
Poisson brackets of Af with the total volume function, 40 ) = 2 i A i(x), J d 3 y \J det E(y)}. 
(This identity holds only if no modulus signs are used around the det A-term, otherwise one 
has to include a factor of sign(det E) in the equation.) In the quantization one substitutes 
the Poisson brackets by commutators. Again the volume operator plays a pivotal role. 

Another role of the operator det E was pointed out recently [17]. This is related to the 
fact that one has to impose a constraint det E > 0 classically, if the basic variables are chosen 
to be the Yang-Mills conjugate pairs (A l a , Ef). This constraint distinguishes the phase space 
of gravity already at the kinematical level from that of a gauge theory. In the lattice theory, 
one can show that all non-vanishing eigenvalues of det FI come in pairs of opposite sign [17], 
and a quantum analogue of det E > 0 can be imposed. In practice, this requires an explicit 
diagonalization of the operator det A 1 . 

Finally, another application is the possible inclusion of a cosmological constant term of 
the form A J d 3 x det E{x) in the Hamiltonian, with the integral taken over the entire spatial 
manifold. There are suggestions that this may be necessary in order to construct a non-trivial 
continuum limit of the (discretized) theory [18]. 

Our present knowledge of the spectrum of det E or the volume operator is only partial. 
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It was observed in [9,10] that its spectrum is discrete and non-zero contributions can only 
come from intersections of the Wilson loop states it acts on. Its diagonalization can in a so- 
called spin-network basis [19] be reduced to a diagonalization on finite-dimensional subspaces 
of Hilbert space. In [7], we gave a general proof for why intersections have to be at least 
four-valent in order to give a non-trivial contribution. Part of the non-vanishing spectrum for 
such four-valent intersections was first given in [8]. These calculations were confirmed in [11], 
where also a general formula for the matrix elements of the volume operator was derived. 
These were given in terms of an orthogonal basis, obtained by decomposing spin-network 
states (closely related to Wilson loop states) with n-valent intersections into three-valent 
ones. Another general expression for matrix elements, with respect to a different orthogonal 
basis, was given in [14], together with formulae for spectra for simple special cases of classes of 
four-valent intersections. (In order to avoid confusion, it should be pointed out that different 
volume operators may differ by overall factors and modulus signs (see, for example, the 
comments in [13]). Still, results on their spectra tend to be closely related.) 

As we mentioned above, det E can be diagonalized on finite-dimensional subspaces. These 
are given by fixing the flux line (or spin) assignments on the edges incoming at a given intersec¬ 
tion where det E acts. Another way of saying this is that there exist operators associated with 
these edges that commute with det E. On the lattice, these are of the form Ei(l)E ai (l) 
(no sum over a) with spectrum yfj{j + 2)/4, j = 0,1,2,..., and related to the measurement 
of area in the a-direction [20]. 

One does not expect to be able to derive analytic formulae for the spectrum of det E as 
the dimension of the finite-dimensional subspaces obtained by fixing the quantum numbers 
j for the incoming edges grows. The explicit evaluation of the spectrum is therefore limited 
by one’s computational ability to diagonalize large matrices. A suitable choice of basis in 
which the matrix elements are evaluated will eventually also be determined by the form of 
the Hamiltonian. 

In this paper we construct further local operators, related to the geometry of the intersec¬ 
tion of a spin network state, that commute with det E, and lead to a further reduction of the 
finite-dimensional Hilbert spaces into smaller subspaces on which det E can be diagonalized 
separately. In mathematical terms, we will be analyzing the irreducible representations of 
a discrete local symmetry group of the quantum states, and decompose the Hilbert space 
accordingly. The method is general and can also be applied to other operators, most im¬ 
portantly the Hamiltonian constraint. This leads to a considerable simplification, as we will 
demonstrate for a class of six-valent intersections, with some of the resulting “superselection 
sectors” becoming totally non-degenerate with respect to their volume spectrum. 
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Our analysis will take place on a cubic lattice, and is therefore valid for all quantum 
states that can be realized on connected subsets of edges from such a lattice. The discrete 
symmetry group is in this case the 24-element group O, the so-called cubic or orthogonal 
group [21]. The analysis may be generalized to intersections with a different geometry and 
symmetry group. 

In the next section, we introduce a labelling of local spin network states (more precisely, 
a labelling for the so-called intertwiners at an intersection) with a simple transformation 
behaviour under the cubic group. For the construction of the superselection sectors, the 
relevant discrete group is a subgroup O ^ of the cubic group times the Z 2 -factor associated 
with total spatial reflection. We describe its irreducible representations and give an explicit 
construction for states transforming according to a given irreducible representation from the 
elements of O ^ x Z 2 -orbits. In Section 3 we give the explicit form of the operator det E 
associated with a vertex (or intersection) n. For a special class of six-valent intersections, 
corresponding to flux line assignments j = (j. j, j, j, j. j ), we perform the orbit decomposition, 
the construction of the superselection sectors, and diagonalization of detFi up to j = 10. All 
eigenvalues, together with their multiplicities in the individual sectors, are listed in Table 9. 
We end in Section 4 with a summary of our results and a brief discussion of the condition 
det E > 0. 


2 Action of the orthogonal group on the quantum state space 

Our kinematical quantum state space is the gauge-invariant sector of the Hilbert space of 
an SU (2)-lattice gauge theory in the Hamiltonian formulation [22]. Following the philosophy 
of spin network states [19], we use flux line labels j = 0,1, 2,... for lattice edges or links l. 
As an overcomplete set of labels for the contractors or intertwiners situated at a vertex n of 
the lattice we use 9 numbers jik, i, k = 1,2, 3 [17]. The index i denotes the incoming edges at 
n from the 1-, 2-, or 3-direction (with respect to some fixed orientation for the three lattice 
axes), and the index k the corresponding outgoing edges. That is, ji 2 denotes the number 
of (unoriented) flux lines routed through n between the (-1)- and the 2-direction, etc. It is 
convenient to arrange the set of nine numbers into matrix form, 


/ in i 12 ii 3 \ 

J '■= ( i21 J22 423 | 

V j.31 J32 J33 J 


( 2 . 1 ) 


Given flux line assignments ji, i = —1, —2, —3,1,2,3, for the in- and outgoing links, it 
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is easy to generate all allowed intertwiner configurations J. The elements of the rows and 
columns of J simply have to add up to the appropriate ji, for example, Y^=i3i,i = i-ij 
j t ,\ = j \. General elements of the cubic group O do not map intertwiners of the form 
(2.1) into themselves (but lead to configurations that by virtue of Mandelstam identities 
can be re-expressed as linear combinations of J-configurations). However, a six-dimensional 
subgroup Cb 6 ) of O leaves the label set invariant. Apart from the identity transformation, 
the non-trivial elements of this subgroup Cb 6 ) act on the J’s according to 


Ri(J) ■■= 


331 

3 21 \ 

( 3 33 

3 23 

3 13 \ 

322 

3 12 

3 32 

333 

323 J 

, R 2 (J) := 1 j 3 2 

322 

3 12 1 

, Rs(J) ■= I J21 

in 

J31 

332 

322 / 

V J31 

321 

in/ 

\323 

J'l3 

333 


Si(J) := 


322 3 23 321 

332 333 331 

J 12 j 13 ill. 


/ 3 33 3 31 3 32 

S2(J) '■= I 313 311 312 

\ 323 321 322 . 


( 2 . 2 ) 



11 

Ri 

R‘2 

Rs 

Si 

S 2 

1 

1 

R\ 

R‘2 

Rs 

Si 

S 2 

R\ 

R\ 

11 

Si 

s 2 

r 2 

Rs 

R-2 

R.2 

S 2 

11 

Si 

Rs 

Ri 

r 3 

Rs 

Si 

S 2 

11 

Ri 

R'2 

Si 

Si 

Rs 

Ri 

R-2 

s 2 

11 

s 2 

s 2 

R2 

Rs 

Ri 

11 

Si 


Table 1 Multiplication table for the subgroup O ^ of the octagonal group. 


The multiplication table for the group Cb 6 ) is given in Table 1. We will also need the total 
space reflection T, corresponding to the transform of the matrix J, 


3 li 

J21 

331 \ 


T(J) := ii2 

322 

332 I • 

(2.3) 

Vii3 

323 

333 ) 



5 




Since T commutes with all elements of (D^\ adjoining it we obtain a 12-element group O ^ x 
T = 0 ( V x 7L-i. This group is important because it is a subgroup of the classical invariance 
group of the lattice function (detE) 2 (see Section 3 below). There is still a redundancy in 
the set of allowed J’s which is associated with Mandelstam constraints and consists of all 
identities of the form 


/ in + i ii2 

( J'21 J22 + 1 

V 1.31 J32 

in ii 2 +1 

321 + 1 J22 

321 J32 

in ii 2 

+ I i21 + 1 i22 

V J31 322 + 1 


il3 \ fj 11 + 1 
J23 1 - I hi 

133 + 1 / V 331 

il3 \ / ill 

J23 j + I hi 
133 + 1 / V i31 + 1 

ii3 + i \ / iii 

J23 J - I J21 
J33 / V hi + 1 


ii2 ii3 \ 

322 322, + 1 | 

32,2 + 1 133 / 

il2 + 1 il3 \ 

322 322 + 1 I 

322 322 J 

il2 il3 + 1 \ 

122 + 1 J23 J = 0. 

322 322 / 


(2.4) 


As already mentioned in the introduction, our aim is to identify the irreducible repre¬ 
sentations of the discrete group O ^ ( O ^ x TL-i), and construct corresponding superselection 
sectors on which the operator detE ((detE) 2 ) can be diagonalized separately. The group 
O^ contains three conjugacy classes of elements namely, {11}, {Ei, R 2 , R 3 } and {Si,S2}. 



{H} 

{Ri} 

{Si} 

{T} 

{TRi} 

{TSi} 

At 

1 

1 

1 

1 

1 

1 

At 

1 

-1 

1 

1 

-1 

1 

E+ 

2 

0 

-1 

2 

0 

-1 


1 

1 

1 

-1 

-1 

-1 


1 

-1 

1 

-1 

1 

-1 

E~ 

2 

0 

-1 

-2 

0 

1 


Table 2 Character table for the group O ^ x Z 2 . 


Following [21], one establishes the existence of three irreducible representations: two one- 
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dimensional ones (called A\ and A2) and one two-dimensional one (called E). They can be 
identified by the values of their characters, i.e. the traces of the matrices representing the 
group elements (which only depend on the conjugacy class). The enlarged group O ^ x T 
has six conjugacy classes and six irreducible representations, since each of the previous re¬ 
presentations gives rise to one of positive and one of negative parity, denoted by a subscript 
+ or —. The characters for the group O ^ x T are given in Table 2. 

In order to establish the contents of irreducible representations of some general represen¬ 
tation, one can make use of the following character formula. It relates the multiplicity m R 
of a given irreducible representation R in a general representation M to the number n R of 
elements in the conjugacy class K, the traces Xic °f matrices in the representation A4 in 
class K , and the characters Xk of the irreducible representation R , 


1 

m R = - 
a 


u kXk 1 Xk- 


(2.5) 


For the group 0^\ d = 6, and for the group x Z 2 , d = 12. 

The possible orbit sizes that occur under the action of the group O ^ x Z 2 on states 
of the form (2.1) are 1, 2, 3, 6 and 12. We distinguish between parity-even and parity-odd 
orbits. In the former, all elements are parity-invariant, i.e. T(J) = J, whereas in the latter 
T(J) ^ J for all J. The contents of irreducible representations of the individual orbit types 
is listed in Table 3 (the numbers in brackets denote the number of J-states). 


1-d even 

^(1) 

2-d odd 

A+(l), Af(l) 

3-d even 

A+(1),E+ (2) 

6-d odd 

A+(l), A~(1),E+ (2 ),E~ (2) 

6-d even 

A+(1), A+(1), E+ (4) 

12-d odd 

At (1), A~ (1), A+ (1), A~ (1), E+ (4), E~ (4) 


Table 3 Representation contents of the O ^ x Z 2 -orbits. 
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Instead of diagonalizing a maximal set of commuting operators built from the elements 
of O (6) x Z 2 (for example, given by {1, R\ + R 2 + R 3 , Si + S 2 , T}) in a given finite-dimensional 
sector of Hilbert space, we have found it computationally simpler to first construct the orbits 
and from those the states transforming according to a definite irreducible representation. 
One way of how this may be done is illustrated in Table 4 (by J we denote in this context 
an arbitrary, fixed element of an orbit). Note that Restates change sign whenever one of 
the Ri is applied. The prescriptions for the instates in the two-dimensional representations 
are of course non-unique. Remember also that (2.4) induces a residual redundancy in the 
sets of states constructed according to Table 4, leading to relations among elements from 
different orbits. For the explicit calculations of the next section (which were performed using 
Mathematica on a DEC AlphaStation 255 4/232 with 64 Mb RAM) these could be handled 
without particular problems. 



all orbits: (1 + T)(J + Ri J + R 2 J + R 3 J + SiJ + S 2 J ) 


all odd orbits: (1 - T )(J + R\J + R 2 J + R 3 J + S X J + S 2 J) 

At 

6-d even orbit: J + S\J + S 2 J — R\J — R\S\J — R\S 2 J 

12-d odd orbit: (1 + T)(J + Si J + S 2 J - R\J - R&J - R\S 2 J) 

A 2 

12-d odd orbit: (1 + T)(J + 5i J + S 2 J - R\J - R\S\J - R\S 2 J) 

E+ 

3-d even orbit: J — S\J , J — S 2 J, 

6-d odd orbit: (11 + T)(J - Si J), (1 + T)(J - S 2 J ) 

6-d even orbit: (1 - Si)(J + R\J), (1 - So)(J + R\J), (1 - S 1 ){R 2 J + S 2 J), 
(1- S 2 ){R 2 J + S 2 J) 

12-d odd orbit: (1 + T)(t - Si)(J + R\J), (11 + T)(ll - S 2 )(J + R\J), 

(1 + T)(ll - Si)(R 2 J + S 2 J), (11 + T)(l - S 2 )(R 2 J + S 2 J) 

E~ 

6-d odd orbit: (11 - T)(J - Si J), (1 + T){J- S 2 J) 

12-d odd orbit: (1 - T)(t - Si)(J + Ri J), (1 - T)(t - S 2 )(J + i?i J), 

(1 - T)(t - Si)(R 2 J + S 2 J), (11 - T)(l - S 2 ){R 2 J + S 2 J) 


Table 4 How to construct states transforming in a given irreducible representation from the 
elements of O ^ x Z 2 -orbits. 







3 Determining the spectrum of the operator (det E) 

We will now apply the results of the previous section to the discretized version det E(n) 
of the function (1.1), acting at a lattice vertex n. The corresponding self-adjoint quantum 
operator we will call D(n) (for convenience rescaled by a factor of 1 with respect to the 
definition in [8]). In terms of the symmetrized link momenta fh(n, a), it is given as 


D{n) := ^y? 7 ai>& e ljk p i {n,a)p j {n,b)p k in,c), (3.1) 

where 

Pi(n,a) = ^( X l + (n,a ) + X l _(n - 1 a, a)), (3.2) 


and X'±{n,a) denote the left- and right-invariant vector fields on the group manifold associ¬ 
ated with the link l = (n, a), with commutators [X±,X±] = ±e ijfc XjL. 

The key observation is that the classical lattice function D(n) = det E{n) is invariant 
under the action of the 24-elenrent cubic group O [21] (whose elements can be thought of as 
the orientation-preserving permutations of the three (oriented) lattice axes meeting at the 
intersection n). By contrast, D(n) changes sign under the total space reflection T (i.e. under 
simultaneous inversion of the three axes). 

As a result of this classical symmetry, eigenstates of D(n ) can be classified according 
to the irreducible representations of O. This set-up is familiar to lattice gauge theorists, 
because it has been employed in analyzing the glueball spectrum of the Hamiltonian in 
four-dimensional SU (3)-lattice gauge theory [23]. As explained in Section 2, in the present 
5C/(2)-context it is convenient to work with a set of local states that is partially gauge-fixed 
with respect to the 0-action, leaving us with O ^ as the residual symmetry group. 

The action of the operator D(n) on states of the form (2.1) can be obtained by considering 
all possible ways in which the triple derivatives contained in Z)(n) can act on sets of flux lines 
routed through the vertex n. Its explicit form is given by 


D 


.711 Jl2 3 13 \ 

.721 J22 J23 | = 

V .7.31 J32 .733 / 


9 



:7ll722733 


"t - T7T j 12 j 2 3 j 31 

lfa 


~T7, j 13 721 7.3 2 
lb 


'(j 11-1 

il2 +1 

il3 \ 

/in - 1 

.712 

.71.3 + 1 \ 

721 

722 - 1 

723 + 1 ) - 

- 1 j ’21 + 1 

322 ~ 1 

723 1 

. V .7.31 + 1 

J 32 

733 - 1 / 

V .7.31 

732 + 1 

733 — 1 / 


in 

712 - 1 

713 + 1 \ . 

f in +1 

ii2 -1 

721 + 1 

722 

723 - 1 ] - | 

721 

722 +1 

731 - 1 

732 + 1 

733 / 

V.731 - 1 

332 

7ii + 1 

712 

jl 3~ l \ 

( in 

712 +1 

721 - 1 

722 + 1 

723 - | 

721 - 1 

722 

731 

732 - 1 

733 + 1 / 

V731 +1 

732 - 1 


+ T^jllj32(jl2 + 2j33 - J13 ~ 721 + 2,722 + 731 + 

lfa 


+ 77jllj23(jl2 - 2j 2 2 - J13 - j‘2 1 - 2j 33 + j 3 l - 

lb 


z 

+ 7^722713(723 + 2j 11 - J21 “ 732 + 2j 33 + jl2 + 
lb 


+ 777.722 J31 (.72.3 
lb 


+ 77733721 (731 
lb 


+ TTlj'S’ij 12 (.7.31 

lb 


- 2j 33 - J21 - 7,32 ~ 2jn + jl2 - 2) 

+ 2j22 - 732 - J13 + 2jl 1 + ,72.3 + 2) 

“ 2jll - J 3 2 - J 1 3 - 2j22 + .723 ~ 


77723 O'll Jl2 + j 12,722 + 712733 + Jl 2 Jl 3 + 2j 12 ) 
lb 


- T7 731 (722 723 + 723733 + jllj23 + 721723 + 2j23) 

lb 


— 77712 (731733 + j 11.731 + 722731 + 731732 + 2j 3 i) 
lb 


+ 77713(711721 + .721.722 + 721733 + 721723 + 2j2l) 

lb 


+ 77721 (722732 + .732.733 + 711732 + .731.732 + 2j 32 ) 
lb 


+ 777.32 (i 1.37.3.3 + .711.713 + 713722 + jl2jl3 + 2jl3) 
lb 
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+ gjllj22(j23 - 

i . . , . 

+ gJllj33Ul2 - 

i . . . . 

+ gj22j33Ul2 — 





/ 

(hi ~ 1 

712 + 1 713 

Vo. 

CO 

to 

+ 

731 - 

-713) | 

, 721 + 1 

722 - 

1 723 





\ 731 

732 

733 




i 

/ jn ~ 1 

712 

713 + 1 

721 

+ 

723 “ 

- 732) | 

72 1 

722 

723 





\ 731 + 1 

732 

733 - 1 




1 

( 711 712 

713 

721 

+ 

731 “ 

-713) | 

, 721 722 

- 1 

723 + 1 





V 731 732 

+ 1 

733 - 1 


+ Tx( Jll(jl2j31 “ 712732 + j 13.723 — i 13.721 + 721723 “ 731732 + .723 ~ 732) + 

16 

.722 (. 721.731 - . 712,713 + , 712.723 - . 721.732 + . 731.732 — , 713.723 + .731 - . 713 ) + 


. 733 ( 712.713 - .721723 + 723731 - 721731 + 712732 ~ 713732 + 712 ~ 721) + 


712713723 - 712713732 + 721723731 “ 713721723 + 712731732 ~ 721731732 + 


/ill 712 713 \ 

712723 + 712731 “ 713721 — 713732 + 723731 — 721732 ) ( 721 722 723 | 

\731 732 733/ 


(3-3) 


Note that the J-states are not orthogonal with respect to the Haar measure on the lattice. 
Still, we will see that within the superselection sectors of the cubic group, degeneracy of the 
eigenvalues of D(n) is largely lifted, which makes orthogonality largely automatic. One can 
verify the following conjugation relations by using the explicit form (3.3) for the operator D: 


RiDRi = D, i = 1,2, 3, SiDSi 


D, i = 1,2, TDT = —D, 


(3-4) 


whence it follows that D obeys the (anti-) commutation relations 


[D, Ri} = 0, i = 1,2, 3, [D, Si + S 2 ] = 0, [D, T\+ = 0. (3.5) 

We conclude that D does not alter the 0++quantum numbers, but maps positive- into 
negative-parity states and vice versa. This latter property suggests a different approach to the 
diagonalization of D(n). Firstly, it follows from [D,T] + = 0 that [D 2 ,T] = 0. Secondly, we 
have shown in [17] that eigenstates of D can be obtained from those of D 2 . For any eigenstate 
X of D 2 with eigenvalue v ^ 0, D 2 (n)x = v 2 Xi one nray construct a pair of eigenstates of 
D(n) with eigenvalues ±w, namely, 

D(x ± T-jDx) = ±WT^ b 2 x + Dx = ±|u|(x ± T~\Dx). (3.6) 

\v u u 
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Since in our search for eigenstates we would like to reduce the size of matrices that have 
to be diagonalized, it is simpler to diagonalize D(n ) 2 first, and then construct eigenstates of 
D(n) using (3.6). The explicit form for D(n) 2 contains 142 terms and can be obtained from 
(3.3). 


One further observation turns out to be useful. Since parity-odd wave functions are 
constructed by weighted sums (with factors ±1) of spin network states, which may sometimes 
vanish, there are always fewer states transforming according to the representations A~, E ~, 
than those transforming according to Af, E + . The most efficient way of diagonalizing D 
(and the one which we follow below) is therefore to start from the wave functions transfor¬ 
ming in one of the negative-parity irreducible representations, diagonalize D 2 , construct the 
images under D of the resulting states (which all have positive parity), and then form linear 
combinations according to (3.6). The number of zero-volume states is given by the difference 
of positive- and negative-parity states. 


j 

■ 

2 

3 

4 

5 

6 

7 

8 

9 

10 

1-d even 

1 

2 

2 

3 

3 

4 

4 

5 

5 

6 

2-d odd 

1 

2 

4 

6 

9 

12 

16 

20 

25 

30 

3-d even 

1 

3 

5 

9 

12 

18 

22 

30 

35 

45 

6-d odd 

0 

1 

4 

9 

18 

30 

48 

70 

100 

135 

6-d even 

0 

0 

1 

2 

5 

8 

14 

20 

30 

40 

12-d odd 

0 

0 

0 

1 

3 

8 

16 

30 

50 

80 

total 

3 

8 

16 

30 

50 

80 

120 

175 

245 

336 


Table 5 Numbers of x Z 2 -orbits for states with j = (j,j,j,j,j,j)- 


As an application of the scheme outlined above, we have analyzed the irreducible represen¬ 
tation contents of the sub-Hilbert spaces with flux line numbers j = (j-i , .7— 2 , J— .3 , Ji , j' 2 , 7.3 ) = 
(j,j,j,j,j,j) (i.e. for genuine six-valent intersections) up to j = 10, and obtained the spec¬ 
trum of D(n). This class of configurations is special in the sense that states with fixed j are 
mapped into themselves by the action of O^ x TLi- The numbers of orbits of a given type 
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for fixed j (before imposing the constraints (2.4)) are listed in Table 5. 


We then proceeded as described above by diagonalizing D(n ) 2 separately on the super¬ 
selection sectors corresponding to different quantum numbers for the action of O ^ x Zo. 
The resulting numbers of linearly independent eigenstates of D(n) with strictly positive and 
vanishing eigenvalues are listed in Tables 6 and 7. 


j 

■ 

2 

3 

4 

5 

6 

7 

8 



A\ 

1 

2 

5 

8 

14 

20 

30 

40 

55 70 


A-2 

0 

0 

0 

1 

2 

5 

8 

14 

20 30 


E 

0 

2 

6 

14 

26 

44 

68 

100 

140 190 


total 

1 

4 

11 

23 

42 

69 

106 

154 

215 290 


Table 6 

Numbers of linearly independent eigenstates of D{n ) 

with eigenvalue > 0. 


j 

■ 

2 

3 

4 

5 

6 

7 

8 

9 10 


total 

3 

7 

12 

19 

27 

37 

48 

61 75 91 


Table 7 

Numbers of linearly independent eigenstates of D{n ) 

with eigenvalue = 0. 


j 

■ 

2 

3 

4 

5 

6 

7 

8 

9 


total 

5 

15 

34 

65 

111 

175 

260 

369 

505 

671 

Table 8 

Number of linearly independent states with fixed flux lines j - 

= {j, 3,3,3,3,3)- 
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Comparing the rows for the A\~, A^- and -E-sectors in Table 6 to the total numbers of states 
in the Hilbert spaces before the O ^ x T-action is taken into account (Table 8), one observes 
that the matrix sizes are reduced considerably when the superselection sectors are treated 
separately. Another interesting feature is that the number of zero-eigenvalues grows less 
rapidly with j than that of the non-vanishing ones (by curve fits we found the dependences 
| j 2 + |j + f for j odd and | j 2 + §j + 1 for j even, as opposed to + § j 2 + 2j + 1). 


Finally, here are our results for the non-negative eigenvalues of the operator D(n) = 
det E(n). The numbers in Table 9 are given in length units, i.e. we have taken the sixth root 
of the eigenvalues of D(n). The three numbers in brackets give the degeneracy of eigenvalues 
in the A \-, A 2 - and E -sectors respectively. 


j = 1 0.820977 (1,0,0), 0 (3) 

j = 2 1.07707 (1,0,0), 0.933429 (1,0,2), 0 (7) 

j = 3 1.28653 (1,0,0), 1.17666 (1,0,2), 1.03890 (1,0,2), 1.02292 (1,0,2), 

0.915413 (1,0,0), 0 (12) 

j = 4 1.46726 (1,0,0), 1.37521 (1,0,2), 1.27163 (1,0,2), 1.25951 (1,0,2), 

1.16697 (1,0,0), 1.12099 (1,1,4), 1.09850 (1,0,2), 1.01217 (1,0,2), 0 (19) 

j = 5 1.62835 (1,0,0), 1.54767 (1,0,2), 1.46030 (1,0,2), 1.45180 (1,0,2), 

1.37044 (1,0,0), 1.35016 (1,1,4), 1.33173 (1,0,2), 1.25452 (1,0,2), 

1.19678 (1,0,2), 1.19140 (1,1,4), 1.16425 (1,0,2), 1.10871 (1,0,2), 

1.09057 (1,0,2), 0.978018 (1,0,0), 0 (27) 

j = 6 1.77500 (1,0,0), 1.70234 (1,0,2), 1.62518 (1,0,2), 1.61898 (1,0,2), 

1.54569 (1,0,0), 1.53375 (1,1,4), 1.52004 (1,0,2), 1.44981 (1,0,2), 

1.42448 (1,0,2), 1.41907 (1,1,4), 1.39619 (1,0,2), 1.34168 (1,0,2), 

1.32815 (1,0,2), 1.26214 (1,1,4), 1.25364 (1,1,4), 1.23241 (1,0,0), 

1.22270 (1,0,2), 1.18305 (1,1,4), 1.15815 (1,0,2), 1.06919 (1,0,2), 0 (37) 

j = 7 1.91050 (1,0,0), 1.84387 (1,0,2), 1.77395 (1,0,2), 1.76921 (1,0,2), 

1.70204 (1,0,0), 1.69398 (1,1,4), 1.68357 (1,0,2), 1.61870 (1,0,2), 

1.60414 (1,0,2), 1.59947 (1,1,4), 1.58195 (1,0,2), 1.52832 (1,0,2), 

1.51891 (1,0,2), 1.48980 (1,1,4), 1.48107 (1,1,4), 1.45463 (1,0,2), 

1.43533 (1,0,0), 1.41349 (1,1,4), 1.39331 (1,0,2), 1.32314 (1,0,2), 

1.32056 (1,1,4), 1.31437 (1,0,2), 1.30965 (1,1,4), 1.27555 (1,0,2), 

1.25331 (1,0,2), 1.24732 (1,1,4), 1.21785 (1,0,2), 1.16207 (1,0,2), 

1.14257 (1,0,2), 1.02552 (1,0,0), 0 (48) 

2.03705 (1,0,0), 1.97517 (1,0,2), 1.91074 (1,0,2), 1.90697 (1,0,2), 

1.84463 (1,0,0), 1.83874 (1,1,4), 1.83056 (1,0,2), 1.77001 (1,0,2), 

1.76040 (1,0,2), 1.75651 (1,1,4), 1.74290 (1,0,2), 1.69067 (1,0,2), 
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j = 8 




1.68385 (1,0,2), 1.66720 (1,1,4), 1.65942 (1,1,4), 1.63881 (1,0,2), 

1.60893 (1,0,0), 1.59606 (1,1,4), 1.58113 (1,0,2), 1.55187 (1,0,2), 

1.54891 (1,1,4), 1.53771 (1,1,4), 1.50965 (1,0,2), 1.50826 (1,0,2), 

1.48282 (1,0,2), 1.47687 (1,1,4), 1.45221 (1,0,2), 1.39757 (1,0,2), 

1.38309 (1,0,2), 1.37801 (1,1,4), 1.37367 (1,1,4), 1.36072 (1,1,4), 

1.32393 (1,0,2), 1.31397 (1,1,4), 1.30466 (1,1,4), 1.28429 (1,0,0), 

1.27155 (1,0,2), 1.23247 (1,1,4), 1.20595 (1,0,2), 1.11432 (1,0,2), 0 (61) 

j = 9 2.15620 (1,0,0), 2.09818 (1,0,2), 2.03812 (1,0,2), 2.03504 (1,0,2), 

1.97663 (1,0,0), 1.97209 (1,1,4), 1.96545 (1,0,2), 1.90848 (1,0,2), 

1.90159 (1,0,2), 1.89833 (1,1,4), 1.88744 (1,0,2), 1.83691 (1,0,2), 

1.83172 (1,0,2), 1.82081 (1,1,4), 1.81415 (1,1,4), 1.79792 (1,0,2), 

1.76329 (1,0,0), 1.75461 (1,1,4), 1.74332 (1,0,2), 1.72777 (1,0,2), 

1.72491 (1,1,4), 1.71478 (1,1,4), 1.69154 (1,0,2), 1.67749 (1,0,2), 

1.66205 (1,0,2), 1.65696 (1,1,4), 1.63812 (1,0,2), 1.60819 (1,1,4), 

1.60323 (1,1,4), 1.58998 (1,1,4), 1.58435 (1,0,2), 1.57430 (1,0,2), 

1.55793 (1,0,2), 1.54379 (1,1,4), 1.53433 (1,1,4), 1.50616 (1,0,2), 

1.48851 (1,0,0), 1.46554 (1,1,4), 1.44411 (1,0,2), 1.42977 (1,0,2), 

1.42829 (1,1,4), 1.42251 (1,1,4), 1.40779 (1,1,4), 1.37128 (1,0,2), 

1.36866 (1,0,2), 1.36845 (1,1,4), 1.36330 (1,0,2), 1.35667 (1,1,4), 

1.32054 (1,0,2), 1.29967 (1,0,2), 1.29326 (1,1,4), 1.26211 (1,0,2), 

1.20555 (1,0,2), 1.18508 (1,0,2), 1.06412 (1,0,0), 0 (75) 

j = 10 2.26912 (1,0,0), 2.21432 (1,0,2), 2.15784 (1,0,2), 2.15526 (1,0,2), 

2.10012 (1,0,0), 2.09649 (1,1,4), 2.09097 (1,0,2), 2.03703 (1,0,2), 

2.03180 (1,0,2), 2.02905 (1,1,4), 2.02009 (1,0,2), 1.97130 (1,0,2), 

1.96720 (1,0,2), 1.95944 (1,1,4), 1.95376 (1,1,4), 1.94061 (1,0,2), 

1.90385 (1,0,0), 1.89752 (1,1,4), 1.88867 (1,0,2), 1.87925 (1,0,2), 

1.87664 (1,1,4), 1.86783 (1,1,4), 1.84934 (1,0,2), 1.82734 (1,0,2), 

1.81716 (1,0,2), 1.81294 (1,1,4), 1.79834 (1,0,2), 1.78328 (1,1,4), 

1.77840 (1,1,4), 1.76633 (1,1,4), 1.74596 (1,0,2), 1.74079 (1,0,2), 

1.73868 (1,0,2), 1.72121 (1,1,4), 1.71284 (1,1,4), 1.69092 (1,0,2), 

1.66208 (1,0,0), 1.66200 (1,0,2), 1.66017 (1,1,4), 1.65364 (1,1,4), 

1.64853 (1,1,4), 1.63864 (1,1,4), 1.63269 (1,0,2), 1.60428 (1,0,2), 

1.60238 (1,0,2), 1.59918 (1,1,4), 1.58717 (1,1,4), 1.55981 (1,0,2), 

1.55607 (1,0,2), 1.53183 (1,0,2), 1.52550 (1,1,4), 1.49950 (1,0,2), 

1.47744 (1,1,4), 1.47486 (1,1,4), 1.46780 (1,1,4), 1.45154 (1,1,4), 

1.44430 (1,0,2), 1.42912 (1,0,2), 1.42294 (1,1,4), 1.41825 (1,1,4), 

1.41035 (1,0,2), 1.40442 (1,1,4), 1.36572 (1,0,2), 1.35747 (1,1,4), 

1.34758 (1,1,4), 1.32752 (1,0,0), 1.31283 (1,0,2), 1.27374 (1,1,4), 

1.24600 (1,0,2), 1.15189 (1,0,2), 0 (91) 


Table 9 Values for x&, where x is a non-negative eigenvalue of D(n). 


The spectrum is rather complex, and becomes more and more spread out with increasing 
j. There seem to be only three different degeneracy patterns, (1,0,0), (1,0,2) and (1,1,4). 


15 



Every eigenvalue that occurs is already contained in the -invariant (i.e. the A\~) sector, 
and is non-degenerate. It would be interesting to see whether this is also the case for more 
general flux line configurations. Note also that the highest eigenvalue for fixed j is always 
non-degenerate. The “volume deficit” observed in [24] (the fact that all volume eigenvalues 
are systematically smaller than expected from the (Laplacian) area eigenvalues, compared to 
the relation one would obtain for a Euclidean reference metric), persists for higher j. although 
it becomes less pronounced. 


4 Summary and outlook 

We have explained in this paper how the diagonalization of the volume operator in the 
loop representation of quantum gravity can be simplified by taking into account discrete 
symmetries at the intersections of the Wilson loop states. Applied to the case of up to six- 
valent intersections on a cubic lattice, this requires the decomposition of Hilbert space into 
the irreducible representations of a six-dimensional subgroup O ^ of the cubic group. One 
finds a set of local operators that commute with det E, and therefore can be diagonalized 
simultaneously. A further simplification arises when one includes the total spatial reflection, 
and diagonalizes the operator (detE) 2 , from which the eigenstates of detE can be obtained 
using the results in [17]. To demonstrate the viability of the method, we have calculated the 
spectrum for a class of six-valent intersections with flux line numbers (j. j , j. j , j, j), for j < 10. 
We found that all eigenvalues are already contained in the 0^-invariant sector, without 
degeneracy. It is conceivable that this sector is also distinguished on physical grounds, but 
this depends on how the continuum limit of the lattice theory will be taken. 

Let us close with some comments on the condition det E > 0 which, as we have argued 
in [17], should be imposed on the quantum state space. It remains to be seen at which 
stage of the quantization it is imposed most conveniently. Ignoring for the moment the zero- 
eigenvalue states, this condition reduces the dimensionality of the Hilbert space by a factor 
2 X , where x is the number of intersections of the lattice. In addition, one obtains a condition 
on physical operators V, namely, that they should not map out of the subspace of states 
with det E-eigenvalues > 0 (or > 0). A sufficient condition is given by the vanishing of the 
commutator 


[P, sign(det E)\ = 0. 


(4.1) 


However, we have not found a simple form of this condition which would not require the 


16 




explicit knowledge of the eigenstates of det E. Note that (4.1) is a rather strong condition 
which, for example, is not satisfied by the lattice analogues of the area operators defined 
in [25]. A less stringent condition is to require that (4.1) be satisfied in the limit as the 
lattice spacing is taken to zero. This is also sensible from a physical point of view, since the 
condition should be independent of the regularization (for example, the version of the area 
operator corresponding to a pure Laplacian [10,20] obviously fulfills (4.1), but coincides with 
other discretized forms of the area to lowest order in the lattice spacing). 

We do not know whether there is a way to formulate a condition like det E > 0 in the 
continuum theory. The regularization used for the volume operator in [14] forces one to use 
modulus signs around the operator, in order that the square root \J |det E\ is well-defined, 
something not necessary in the case of the lattice regularization. Nevertheless, we think 
that such a constraint on states, along with operator conditions of the type (4.1), should 
be imposed in the quantum theory - even if this leads to new complications - because they 
describe a property of the gravitational theory. 
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